Seeking better times: atomic clocks in the generalized Tonks-Girardeau regime 
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First we discuss briefly the importance of time and time keeping, explaining the basic function- 
ing of clocks in general and of atomic clocks based on Ramsey interferometry in particular. The 
usefulness of cold atoms is discussed, as well as their limits if Bose-Einstein condensates are used. 
We study as an alternative a different cold-atom regime: the Tonks-Girardeau (TG) gas of tightly 
confined and strongly interacting bosons. The TG gas is reviewed and then generalized for two-level 
atoms. Finally, we explore the combination of Ramsey interferometry and TG gases. 

PACS numbers: 



I. INTRODUCTION 

It is hardly necessary to insist on the practical importance of time since we all live attached to a time machine, the 
wrist watch, which organizes our daily routine. For running the Blaubeuren Conference, a precision of one minute 
is enough, but for plenty of activities fundamental to modern society a much more precise time keeping is necessary: 
banks, electric power companies, telecommunications, or the GPS use atomic clocks. 

A clock is a device that, in a sense, "produces" time by counting stable oscillations, for example of a pendulum. 
Clocks of many different types have been used along history by different civilizations but the Earth rotation has been 
for centuries the master clock. Nevertheless, we know that the rotation speed may be affected by a series of physical 
phenomena, such as tidal friction, which slows down the Earth at a non-negligible pace in the geological scale (a 
poor dinosaur had to rush for a sandwich in stressful days of 23 hours); even the accumulation of snow on the top of 
the mountains in the winters of the Northern hemisphere may affect the rotation speed -by conservation of angular 
momentum- making winter days indeed longer than summer days. If we kept the definition of the second attached 
to the actual Earth's rotation, we should also change the laws of physics from summer to winter, and add proper 
correction terms to account for tidal friction. This sounds very unreasonable. Instead, we could also try to correct 
for these perturbing effects "by hand" , but unfortunately not all the perturbations are easy to predict or understand 
(e.g., those due to magma motion), so corrections using other astronomical periods were explored, but they were 
very cumbersome to handle, and not stable enough anyway. Fortunately nature has a lot of stability to offer in the 
opposite direction, in the world of the small. 

Atomic clocks, in particular, count the oscillations of the field resonant with an atomic transition, most frequently 
a hyperfine transition of caesium-133 (which defines officially the second as the time required for 9192631770 periods 
of the resonant microwave field). An external quartz oscillator is locked by a servo loop to a resonance excitation 
curve between two hyperfine states \g) and |e) so that its frequency is always adjusted to the maximum of the curve, 
and thus to the natural frequency of the selected transition. 

In more detail, see Fig. (TJ caesium atoms are heated in an oven to produce a beam. A first magnet filters out 
the excited atoms so that only ground state atoms enter into the microwave cavity in which the field wavevector is 
perpendicular to the atomic motion, and the frequency, locked to an external quartz oscillator, is very close to the 
transition frequency. This excites some atoms which are selected with a second magnet and later detected. The 
excitation probability compared with previous runs at slightly different frequencies tells us how far or close is the field 
frequency to the transition frequency, and this information is used to modify the excitation frequency so that it stays 
as close as possible to the maximum of the excitation curve, see Fig. [2l The clock includes the appropriate counting 
electronics which we shall not discuss here. A very sharp peak is clearly desirable to minimize the fluctuation of the 
external oscillation frequency around the natural one. This contributes to the stability of the clock and explains why 
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FIG. 1: Principle of an atomic clock based on a "Rabi", or one-cavity scheme. Deflection angles are greatly exaggerated. Actual 
clocks operate usually with two separated cavities ( "Ramsey scheme" ) . 
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FIG. 2: An excitation probability curve versus detuning (arbitrary units). 



a beam configuration is chosen: a perpendicular excitation with respect to the atomic motion avoids the Doppler 
effect and its associated line broadening. 

In actual clocks things are slightly more complicated because instead of a single cavity (Rabi scheme) there are 
actually two (Ramsey scheme). The reason for the two cavities is that they produce a quantum interference between 
two possible paths corresponding to excitation in the first cavity or excitation in the second, and this interference may 
be used to sharpen the resonance and to make it less dependent of inhomogeneities of the fields. The Hamiltonian 
describing the process is 

^ = ^- + \^{m9m + |e)(5|), (1) 

where Vt is the Rabi frequency and J the detuning (laser frequency minus transition frequency). In the semiclassical 
approximation the center of mass motion is classical so that p and x become numbers: p is the momentum and x a 
linear function of time; T is the flight time between the fields and r the time at each field region, see Fig. [3l Note 
that the second path (excitation in the second field) does not pick up any extra phase during free flight between the 
fields but the first one does (e^'^^). The consequence is a final excitation probability proportional to cos^(^T/2), which 
makes clear that longer free- flight times are desirable in order to achieve a narrower central fringe. 

This motivates the use of cold atoms in time-frequency metrology. They are slow and big T's may be produced, but 
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FIG. 3: The two main interfering paths in Ramsey's method. Here |1) = |^), |2) = |e) 



there are other advantages: first they may reduce dramatically velocity broadening (the averaging of the fringes due 
to different velocities and times), and also they may lead to fundamentally new effects with coherent N-body states 
(e.g. entanglement has been proposed to beat quantum projection noise limit). 

So why do not we try to use Bose-Einstein condensates? Before we answer this question, let us tell a little story to 
show that "the colder the better" is not an infallible motto. The Chinese made very complex water clocks centuries 
ago in which the falling water moved a wheel with small cups. Some versions of these clocks were also built in Europe, 
but they obviously had serious problems in cold winters: the water froze and the clock did not work. 

As in our water clock example, the improvements associated with low velocities and narrow velocity distribution in 
a Bose-Einstein condensate may be compensated by negative effects, such as collisional shifts and instabilities leading 
to the separation of the gas cloud [J, Q. So BECs do not seem to make good clocks, but there are other cold-atom 
coherent regimes. In particular, the Tonks- Girardeau (TG) regime of impenetrable, tightly confined Bosons subjected 
to "contact" interactions [Si, i4] is opposite to the BECs in some respects and offers a priori interesting properties for 
metrology: the TG requirement of strong contact interactions, implies similarities between the Bosonic system and 
a "dual" system of freely moving Fermions, since the particle densities of the TG Bosons and the dual system of 
Fermions are actually equal, so we may expect that the strong interaction regime is actually an asset for a good clock. 
Other important feature of the TG gas is its one dimensional (ID) character. Olshanii showed [H, 0, 13 that when a 
Bosonic vapor is confined in a wave guide with tight transverse trapping and temperature so low that the transverse 
vibrational excitation quantum Tiuj^ is larger than available longitudinal zero point and thermal energies, the effective 
dynamics becomes one dimensional, and accurately described by a ID Hamiltonian with delta-function interactions 
giD^{xj — x^), where xj and X£ are ID longitudinal position variables. This is the Lieb-Liniger (LL) model [8J. The 
coupling constant giD can be tuned by varying the magnetic field (and thus the three dimensional s-wave scattering 
length) or the confinement (confinement induced resonances) [5J, 10]) near a Feschbach resonance; it is thus possible 
to reach the Tonks-Girardeau regime of impenetrable Bosons, which corresponds to the giD oo limit of the LL 
model. The first experiments were carried out in 2004 [9, 10], whereas the TG gas model was proposed and solved 
exactly in 1960 by M. Girardeau 0,0] exploiting the similarities between the Bose gas and its dual ideal Fermi system 
mentioned above with the so called Fermi-Bose mapping. In actual confined gases the TG regime requires a large 
ratio between the chemical potential and the kinetic energy, 7 = mgiD/h'^n [il|, n being the linear density. Note that 
a smaller density is favorable (to make this more intuitive, think of a rescaling of the lenghts in which the interatomic 
interaction is effectively closer to a delta function as the gas becomes more dilute). On the other hand we do not 
want n to be too small for a clock since the total number of particles (and thus the signal) would also decrease. 

The recipe to construct an TV-boson TG-wavefunction by Fermi-Bose mapping is as follows: 

• Build a Slater determinant for ideal (free) spinless Fermions using one-particle orbitals 



det^u^^(t)n{xk). 



/TV! 



(2) 



Apply the "antisymmetric unit function" A = Y\i<j<k<N ^d'^i^k — Xj)i as 

V^b(^i, . . ',Xn) = A{xi, . . .,XN)i^F{xi, . . . ,XAr). 



(3) 



In metrology and atomic inter ferometry, the tight ID confinement along a waveguide has pros and cons: the absence 
of transversal excitation and motional branches may lead to an increased signal, but the confinement is by itself 
problematic for frequency standard applications, since the necessary magnetic or optical interactions will in principle 
perturb the internal state levels of the atom. Several schemes have been proposed to mitigate this problem and 
compensate or avoid the shifts [J, [lH, [l^ , and we shall here assume that such a compensation is implemented. 
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To investigate the implications in Ramsey inter ferometry of a strongly interacting ID gas, it is first necessary to 
generalize the TG gas model by adding internal structure. The resulting idealized model remains exactly solvable if 
the collisions produce internal state and momentum exchange. We shall also briefly discuss the Ramsey scheme in 
the time domain, which turns out to give very similar results for reasonable parameters. 



II. TWO-LEVEL TONKS-GIRARDEAU GAS WITH EXCHANGE, CONTACT INTERACTIONS 



The first step to generalize the TG gas to one with internal states is to use instead of one-particle orbit als, one- 
particle "spinors" 



b=g,e 



(4) 



where n = 1,2, 3. ..is a label to distinguish different spinors (it will later on correspond to states prepared as harmonic 
oscillator eigenstates of a longitudinal trap), and b may be g (ground), or e (excited) {g and e do not necessarily 
correspond to states with definite values of the component of the electronic spin in one direction.) Two-particle states 
may also be formed as 



(5) 



b,b' 



and similarly for more particles. In \bb^), b is for particle 1 and b' for particle 2. This will in some equations be 
indicated even more explicitly as 6i, 62, etc. 

Analogously to the steps to construct the ordinary TG gas wavefunction, let us first define a Fermionic state for 
two noninteracting particles with internal structure as a generalized Slater determinant. 
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More explicitly. 
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Weg). (6) 

These Fermions do not interact among themselves, but they could interact with an external potential. 

A Bosonic wavefunction of interacting atoms, symmetric under [xi ,bi) ^ {xj , bj ) permutations may be now obtained 
by means of the Bose-Fermi mapping, '^b{xi,X2) = f{xi, X2), where the antisymmetric unit function is ^ = 
sgn(x2 -xi). 

For the sector xi < X2, 



2^/Hb{xuX2) = 
+ [4'\X,)4'\X2) 

+ [4^\x^)42\x2) 
+ [ct>['\xi)4'\x2)]\ge) 
+ [ct>'f\x^)4'\x2)]\eg) 
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(7) 
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FIG. 4: Collisions for triplet and singlet potentials. 



In the complementary sector xi > X2 we would get the same form except for a global minus sign. The resulting Bosonic 
state is discontinuous at contact so it cannot represent noninter acting Bosons. The generalization to A/'-particles is 
straightforward [20*]. 

To understand the physical meaning of the contact interactions implicit in Eq. ([7j), let us now introduce Pauli 
operators, 

= |^)(e| + |e)(^|, 
ay = z{\g){e\-\e){g\), 

= |e)(e|-|^)(^|, (8) 

and the corresponding 3-component operator Sj = aj/2 for particle j analogous to the spin- 1/2 angular momentum 
operator. For two particles S = S1 + S2, and has eigenvalues S{S-\-l) with S = and S = 1 corresponding to singlet 
and triplet subspaces. These subspaces are spanned by | — ) = {\eg) — \ge))/\^ and {\gg)^ |ee), |+) = {\eg) + \ge))/V^} 
respectively. 

Assume now the Hamiltonian 

t.2 2 

Here X12 = Xi — X2, P12 = | — )(— I = | — 81-82 is the projector onto the subspaces of singlet functions, and 

^12 — \99){99\ + \^^){^^\ + is a projector onto the triplet subspace. 

The internal Hilbert space can be written as the sum of singlet and triplet subspaces as Tig^Ht- Suppose that the 
reflection amplitude for relative motion in such representation takes the values +1,-1 in singlet and triplet subspaces 
respectively. The particles are impenetrable and these values correspond to a hard wall potential vt at X12 = for all 
collisions in triplet channels \gg) \gg)^ |ee) \ee) or |+) |+), whereas Vs has, in addition to the hard wall at 
X12 = 0, a well of width / and depth V, so that the reflection amplitude becomes = +1 in the limit in which the 
well is made infinitely narrow and the weh infinitely deep, keeping (2mVyo/fi^Y^'^L^ = 7r/2 [11, [H, [H, [T^ • (In the 
Fermionic Tonks- Girardeau gas the well applies to the triplet subspace and not to the singlet subspace as here.) 

If we translate the above into the g, e-basis and for the sector xi < X2 this implies that in all collisions between 
atoms with well defined incident momenta, they interchange their momenta (the relative momentum changes sign), 
as well as their internal state, with the outgoing wave function picking up a minus sign because of the hard-core 
reflection, as shown in Fig. [5l For xi < X2 and equal internal states such collision is represented by 

eikx^^ik'x2\^^^ _ [e^f^'^ie'^'^'jlbb), (10) 

whereas for h ^h' ^ 
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FIG. 5: Diagrammatic representation of the collisions for particles with different internal states in the sector xi < X2 (the last 
two terms in Eq. (U}): the particles interchange their momentum and internal state picking up an additional phase (minus 
sign). 
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FIG. 6: Schematic setup for Ramsey interferometry of guided atoms in the spatial domain, 
ground state and the probability of excitation is measured after passing the two fields. 
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For equal internal states the spatial part vanishes at contact, xi = X2, whereas in the non-diagonal case it does not, 
but in Eq. (pTj) only the "external region" is considered, disregarding the infinitely narrow well region. 

For an implementation of these contact interactions, undesired inelastic collisions should be suppressed by con- 
finement at low collision energies [l^, which will also produce effective triplet reflection coefficients close to — 1, 
independently of the internal state. 

For two Bosons in the singlet subspace, the space wavefunction is antisymmetric, so that s-wave scattering is 
forbidden; therefore the interactions are governed to leading order by a 3D p-wave scattering amplitude and can be 
enhanced by a ID odd- wave confinement-induced Feschbach resonance (CIR), which allows in principle to engineer 
Vs and achieve a strong attraction as required above. 



III. THE RAMSEY INTERFEROMETER WITH A TG GAS 



We consider as the initial state for the Ramsey experiment N two- level atoms in the Tonks- Girardeau regime, 
confined in their ground internal states in a harmonic trap of frequency co. The cloud is released by switching off the 
trap along the x-axis at time t = (cj = for t > 0); a momentum kick hko is also applied, so that it moves towards 
the two oscillating fields localized between and / and between / + L and 21 ^ L (Fig. [6|). The central initial position 
of the harmonic trap xq < is such that \xq\ » 6n = [{N + l/2)fi/(mcc;)]^/^, to avoid the overlap with the spatial 
width (root of the variance) of the highest occupied state. 

In an oscillating-field-adapted interaction picture (which does not affect the collisional Hamiltonian) and using the 
Lamb-Dicke, dipole, and rotating- wave approximations the Hamiltonian is given, for each of the particles, by Eq. ([1]), 
For the explicit x dependence of Q{x) we assume Q{x) = ^} for x G [0, /] and x e [/ + L, 2/ + L] and zero elsewhere. 
We include the interparticle interactions implicitly by means of the TG wave function and its boundary conditions at 
contact. 

The Ramsey pattern is defined by the dependence on the detuning of the probability of excited atoms after the 
interaction with the two field regions. For the TG gas it follows that the total excitation probability is obtained as 
the average of the excitation probabilities for each spinor, Pg = ^n=iP^^ - O^ce a particle incident from the left 
and prepared in the state e*^°'^^~^°Vn(^ ~ ^o)|^) at t = has passed completely through both fields, the probability 
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amplitude for it to be in the excited state is 

where 4>n{k) is the wavenumber representation of the kicked n-th harmonic eigenstate, 



<Pn{k) 
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g-5e(fc-feo)%-ifca;o^^[y2,5o(/c - fco)], 



(12) 



(13) 



the momentum in the excited state is g = ^^/c^ + 2mA/fi, the spatial width of the n = state is Jo = [^/ (2ma;)]^/^, i^n 
the Hermite polynomials, and Tg^ is the "double-barrier" transmission amplitude for the excited state corresponding 
to atoms incident in the ground state (the excited state probability for monochromatic incidence in the ground state 
is ||T^eP)- A fully quantum treatment of Tg^ can be found in [l^. We have done numerical simulations for / = 1 cm, 
L = 10 cm, N = 10, and vq = 1 cm/s. The variation of the excitation probability for different harmonic eigenstates is 
negligible, and the curves for the central fringe are indistinguishable from the semiclassical result of Ramsey (classical 
motion for the center of mass, uncoupled from the internal levels). 



n'r\ fAT\ A fQW\ . /AT 

cos — -— sm sm 

2 J \ 2 J ft' V2y \ 2 



(14) 



where r = I/vq, T = L/vq, and = {Q^ + A^Y^^. 

The broadening of the central fringe by increasing n due to the momentum broadening of vibrationally excited 
states is quite small for the few-body states of our calculations, N = 10, which is in fact of the order of current 
experiments with TG gases {N ^ 15,50 in [i,ll3), see [20] for details. 

We should also examine the quantum projection noise due to the fact that only a finite number of measurements are 
made to determine the excitation probability. The error to estimate the atomic frequency from the Ramsey pattern 
depends on the ratio [21j 
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(15) 



which we have calculated at half height of the central interference peak. Even though r is theoretically smaller for the 
TG gas than for independent atoms, for reasonable parameters the ratio r essentially coincides with that for freely 
moving, uncorrelated particles ^] and, for L >> / it gives l/{T^/N) for all 5 [20]. 



IV. RAMSEY INTERFEROMETRY IN THE TIME DOMAIN FOR GUIDED ATOMS 

An alternative to the previous set-up is the separation of the fields in time rather than space but, at variance with 
the usual procedure, keeping the gas confined transversally at all times as required for the ID regime of the TG 
gas. Fig. [71 This configuration is thus similar to the one used in ion trap clocks. Because of the tight confinement 
the transverse vibrational excitation is negligible so that the Ramsey pattern is given by the standard semiclassical 
expression irrespective of the value of n. The whole TG gas therefore produces the usual Ramsey pattern The 
advantage with respect to trapped ions is the augmented signal, and a disadvantage would be the need to take care 
of, and possibly correct for, perturbing effect of the confinement on the levels of the neutral atoms. 



V. SUMMARY AND DISCUSSION 



Our model of a generalized Tonks- Girardeau regime suggests that very cold atoms may provide good atomic clocks 
if strongly confined. The advantages are: a very large time between field interactions, no collisional shifts, no velocity 
broadening, and no instabilities as in BECs. More work is needed to ascertain the role of strongly interacting gases 
in interferometry and time-frequency metrology: to go from models to actual atoms, and to compensate or eliminate 
shifts due to confinement forces. 
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FIG. 7: Schematic setup for Ramsey interferometry in time domain. The TG gas is confined in a cigar-shaped trap and 
illuminated by a laser in ^/-direction. 
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